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We find a specific meclianism of background spin currents in two-dimensional electron systems 
with spatially nonuniform spin-orbit interaction (SOI) at tiiermodynamic equilibrium, in particular, 
in the systems consisting of regions with and without SOI. The spin density flows tangentially to 
the boundary in the normal and SOI regions within a layer of the width determined by the SOI 
strength. The density of these edge currents significantly exceeds the bulk spin current density. The 
spin current is polarized normally to the boundary. The spatial distribution of equilibrium spin 
currents is investigated in detail in the case of a step-like variation of the SOI strength and potential 
at the interface. 



The generation and manipulation of electron spin 
polarization in semiconductor nanostructures solely by 
means of electric fields has attracted significant atten- 
tion as an interesting manifestation of spin-orbit inter- 
action (SOI) as well as a valuable capability for spin- 
tronic devices A key point in solving this problem is 
searching for mechanisms of effective generation of purely 
spin currents. In this connection a great deal of interest 
was paid recently to the dissipationless spin currents in 
semiconductors with SOI'^'^ and particularly, to the spin 
currents in two-dimensional (2D) electron systems un- 
der the thermodynamic equilibrium.—!^ The existence of 
equilibrium background spin currents was first pointed 
out by Rashba^ for an infinite 2D electron gas (2DEG) 
with SOI. This phenomenon has arisen a wide discussion 
of the definition of spin currents^i^ ' '— and the possi- 
bility to observe the equilibrium spin current o^^d^ . The 
definition of the spin currents is known to be somewhat 
arbitrary since the spin is not conserved in the presence 
of SOI. The discussion has led to the evident conclusion 
that observable effects do not depend on any definitions. 
One can use the standard definition of spin current to 
calculate correctly measurable values. Sonini^ proposed 
to detect the equilibrium spin current by measuring the 
mechanical torques near the edges of the SOI medium. 
Sun et alii argued that the equilibrium spin current can 
be interpreted as a persistent current, and considered the 
possibility to determine this current in a ring device by 
measuring an electric field it generates. 

In this paper we draw attention to another aspect of 
the spin current problem. We study equilibrium spin 
currents in a 2D electron system with nonuniform SOI. 
We were initially motivated by a question of whether the 
equilibrium spin currents penetrate from a 2DEG with 
SOI into normal 2DEG. This problem is interesting by 
two reasons: (i) if the spin current could penetrate into a 
normal 2DEG, it could be a measure of the spin current 
in SOI medium since in a normal 2DEG the spin cur- 
rent is unambiguously defined;^"^ (ii) many realistic 2D 
structures usually contain regions in which Rashba SOI 
is induced by a normal electric field produced by gates 
forming mesoscopic structures or by built-in charges. In 



particular, the SOI apparently exists at edges of con- 
strictions in 2DEG and may affect the spin-dependent 
electron transport observed in experimentsi^. 

We study the equilibrium spin currents in a 2D elec- 
tron system consisting of regions with and without SOI to 
come to unexpected conclusion that very effective mech- 
anism of the spin current generation acts near the bound- 
ary. This mechanism generates the equilibrium edge spin 
current of significantly higher density than the bulk spin 
current in an infinite 2DEG with SOI. The edge spin 
current density is proportional to the SOI constant a, 
while the bulk current is proportional to . To be more 
specific, the edge and bulk currents are, respectively, pro- 
portional to akp and o'^, where a is the SOI wave vec- 
tor and kp is the Fermi wave vector. The difference in 
the magnitudes of these currents is caused mainly by the 
fact that the edge spin current is produced by electrons 
in the whole energy range of occupied states, while the 
bulk spin current is generated in an energy layer, the 
width of which equals to the characteristic SOI energy. 
A remarkable result is that the equilibrium edge spin cur- 
rents exist in the normal 2DEG where the spin current 
is strictly defined. 

Let us consider the equilibrium spin currents in a 2D 
system consisting of semi-infinite regions with and with- 
out Rashba SOI (Fig.[T]). The SOI is described by a step 
function a{x) — aQ{—x). The potential energy U{x) is 
also step function, the potential in the SOI region being 
higher than in the normal one. This arrangement of the 
potential shape allows one to study effects of both prop- 
agating and evanescent modes in the SOI region on the 
spin state in the normal region. 

First, we briefly describe the spectrum of electron 
states in the SOI region. There are three types of modes 
studied in detail in Ref. [isl . The eigenfunctions are 

^^^^C^j^XA(k)^^,(,^,+,^^)^ (1) 

where A stands for spin states, k = [kx,ky), and C is 
a constant. Since in the case being considered, ky is a 
conserving value, we will drop this index in the notations 
of wave functions. The characteristic wave vector and 
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FIG. 1: (Color online). The energy diagram of the normal- 
SOI structure. £(|fc|) denote electron spectra in the SOI re- 
gion; (+) and (-) are upper and lower states with propagat- 
ing and decaying evanescent modes; (o.e.s.) is "oscillating" 
evanescent modes. Dotted is the energy layer, in which the 
equilibrium spin currents are generated in the bulk of the SOI 
reservoir. £o(g) is the spectrum in the normal region. 



energy of the SOI are a ~ majY? and Sso = h^a^ /{Im). 
In the energy range e > —Ego, there are two modes: 

(i) propagating modes |A, fc^), with kx being real, ~ 
k'^; A = ±; and XA(k) = \{ky + ik'x)/k; 

(ii) "decaying" evanescent modes |A, fc"), with kx being 
imaginary, k^ = ik'J.; A = ±; and XA(k) = X{ky — k'^)/k. 

The energy of both modes above is 



eA,k 



2m 



xky 



(2) 



where k = {kl+ky Y^"^. At a given energy e > 0, the wave 
vector k takes two values fci and k2 corresponding to the 
states with positive and negative chiraUties (A = ±): 



-Aa -I- \j2melW- + 



(3) 



In what follows one needs to know which states exist 
at a given energy. For e > 0, three cases are possible. If 
\ky \ < ki, there are two pairs of propagating states \zLk[^) 
and I ± k'2x)- li ki < \ky\ < fc2, the states with A = -1-1 
are "decaying" evanescent | ± fc"^.) while the states with 
A = — 1 remain propagating | ± k'2x)- If ^2 < \ky\, both 
states with A = ±1 are "decaying" evanescent | ± fc"^) 
and|±fc^;). 

In the range —Sso < e < 0, the states with positive 
chirality disappear and all four states have A = — 1. In 
this case there are two characteristic wave vectors 



ki,2 = a T \l2melW- 



(4) 



If \ky \ < fci, two pairs of propagating states | ± k'l^) and 
I ± Kx) exist. If fci < \ky\ < k2 



are "decaying" evanescent 

If fc2 < \ky 



± Kx) 



a pair of the states 
and another pair is 
propagating jifej^.). If fc2 < \ky \ all states are "decaying" 
evanescent, | ± fc"^) and | ± fc2a:)- It is worth noting that 
it is just the energy layer that generates the equilibrium 
spin currents in an infinite 2DEG with SOIi^ 



In the energy range e < —ego: the electron states are 
"oscillating" evanescent |fc^,fc"). They are described by 
a complex wave vector kx — k'^+ ik'^, 




(5) 



Here fc^, fc", and ky are coupled by an equation described 
in Ref. 15. At a given energy there are four complex 
vectors kx corresponding to different signs of k'x and fc". 

Note that in the bulk of the SOI reservoir solely the 
propagating states |A, fc^, ky) exist. Evanescent states can 
be realized only near the interface. 

In the normal reservoir the electron states |A, q, ky) are 
propagating. They are degenerate by the spin: 



1 



2U/2 



J(qx+kyy) 



(6) 



where q is the x component of the wave vector. 

Now let us turn to the normal-SOI structure. The wave 
functions can be easily written for a region of a width 
shorter than electron mean free path near the interface. 
The electron states are formed by the waves coming from 
the left and right reservoirs. They partially transmit the 
interface and reflect from it. Correspondingly there are 
four types of wave functions of the form: 
(i) for right moving electrons 



\s,kx)+^r^^s,\s',-kx), x<0, 

s' 

s' 

(ii) for left moving electrons 

s' 



(7) 



(8) 



Here r^^J" and t^^^^ are reflection and transmission ma- 
trices. Index 5 = 1,2 numerates the states in the SOI 
and normal regions, which participate in the scattering 
process at the interface. For the state vectors in the nor- 
mal region, the index s = 1, 2 corresponds to the chirality 
A = ± in Eq. ([6]). For states in the SOI region, the sit- 
uation is more complicated. In the energy region e > 0, 
the index s = 1, 2 also matches with the chirality A = ±. 
In the energy layer —Ego < £ < 0, where A = —1, the 
index s = 1,2 labels two roots defined by Eq. (|4]). For 
£ < —Sso the index s labels two states with negative fc" 
and different signs of fc^: s = 1 corresponds to fc^ > 0, 
and s — 2 does to fci < 0. 
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The reflection and transmission matrices are deter- 
mined by matching the wave functions ([7]) and ([5]) at 
the interface using the standard boundary conditional^: 



= 0, 



dx 



{iaa-y - l3ky(Tz)i^ 



+0 



= 0, (9) 



where (3 takes into account the Rashba SOI caused by an 
in-plane gradient of the potential at the interface. 
The spin current is defined as 



J^{r,s,\^) = ^(V-gl-^^. -t-«.a^|Ck> . (10) 

where r labels the right and left moving states and 7 
stands for the polarization. This definition is strictly vali- 
dated for the normal region of the structure under consid- 
eration, where Eq. (jlOp follows from the spin conservation 
law. In the SOI region the spin current is not conserved 
and therefore can not be unambiguously defined. Equa- 
tion pop represents the commonly accepted and phys- 
ically appealing definition of spin current. ^'ii^i^i^iiSiii 
We emphasize that the spin currents in the normal re- 
gion do not depend on the definition of the spin current 
in the SOI region. 

The total current is a sum over all electron states below 
the chemical potential /i, 
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de dkyD'f\e,ky)j;'{r,s,k), (11) 



where D^s\e, ky) is the density of states for waves inci- 
dent on the interface from the left or right reservoirs. 

Calculations show that in the normal region all com- 
ponents of the spin currents JZ^ directed normally to the 
interface are absent, though in an infinite 2DEG with 
SOI the spin current is nonzero.— Hence, the normal 
component does not come outside the SOI medium, in 
any case through a smooth border line we are consider- 
ing here. The tangential component with tangential po- 
larization is also absent. Only nonzero components 
are and J^y . Thus in the normal region spins fiow 
tangentially to the interface and are polarized normally 
to the interface. The spatial distribution of the spin cur- 
rents is presented in Fig.O In the normal region the spin 
current is located at a distance of the order of 1/a and 
after some oscillations vanishes in the bulk. 

In the SOI region the spatial distribution of the tan- 
gential spin currents essentially depends on the chemi- 
cal potential. If < —£so the spin current is caused 
by the "oscillating" evanescent states and therefore van- 
ishes with distance. In the case ^ > —£so^ the 2;-polarized 
current vanishes with the distance while the x-polarized 
current goes to a finite value coinciding with the spin 
current in infinite systemi^ The normal components 
turn to zero at the interface for all polarizations. Notice 
that the spin density is zero everywhere. 

Of most importance are the following facts: (i) near 
the interface the spin current is much higher than the 
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FIG. 2: (Color online). The distribution of spin current di- 
rected tangentially with spin polarization in (a) x direction 
and (b) z direction. Curves 1, 2, and 3 correspond to the 
chemical potentials fj,i, ^2, and /is shown in Fig. [T] Thin 
horizontal line shows the equilibrium spin current density 
in infinite 2DEG with SOI. The current is normalized to 
= h^a^/^m. Calculations are carried out for U = Gsso 
and /3 = -0.01. 



maximum of the current attainable in an infinite system; 
(ii) the near-boundary spin current is generated by elec- 
trons with any energy in contrast to the case of infinite 
system where the equilibrium spin current is generated in 
the energy layer — e^o < e < 0. These facts evidence that 
near the interface another mechanism of the equilibrium 
spin current acts, which is caused by the spin-dependent 
scattering of electrons by the interface 

We demonstrate this mechanism by considering a sim- 
plified situation where an uniform 2DEG contains a re- 
gion of the SOI in the form of a strip lying along the y 
axis. The SOI strength a is a continuous function of x 
localized at a distance 21. The Schrodingcr equation is 



f^d^ 
2m dx"^ 



Vix)^ 



where 



V{x) = a{x) ( kya^ + iay — 



2m 



i da . 
2dx^ 



(12) 



(13) 



4 



Suppose that the SOI term is small and treat it as a 
perturbation to study the scattering of electrons incident 
from the left and right reservoirs. Unperturbed wave 
functions are 



V^Ajfea; ,ky ^JTj^ ^ Y 



J{k^x+kyy) 



(14) 



where x = {ky + ikx)/k. The scattered wave functions 
are easily found in the Born approximation. We are in- 
terested mostly in the wave functions outside the SOI 
strip. For x > I one finds 



dxa{x) and 



kx <0, 



(15) 



where ag — 

dxa{x) exp{2ikxx). 
Using these wave functions we get simple expressions 
for the spin density and spin currents averaged over the 
equilibrium distribution function. The results are as fol- 
lows. The spin density equals zero, AS = 0. It is curious 
that in each state V'A.k, the spin density is perturbed, but 
summing over A nullifies all spin components excluding 
the X component. The latter is an odd function of ky: 



A5,(k) 



2mky 
h\kx 



-Im 



2i\k^ 



(16) 



This spin density perturbation is produced only by 
backscattered electrons. Transmitting electrons do not 
contribute to the spin density, in the first approxima- 
tion. The integration of AS'a;(k) with respect to k with 
the equilibrium distribution function results in the zero 
AS. 

The tangential spin current in each electron state is 
simply a product of the spin density and tangential ve- 
locity. So the only nonzero component of the spin cur- 
rent, after summing over A, is J^f (k) — {hky/m)ASx(k), 
which is an even function of ky. Finally the total spin 
current is 



V 



2 
3^ 



dk 



k^)^^^ Im [a2ke 



2ikx'\ 



(17) 



where k^ is a low-energy cutoff caused by the violation 
of the Born approximation. Integration leads to result 
qualitatively similar to those shown in Fig. [2ta). The 
z component of the spin current appears in the second 
approximation. 

The mechanism of the equilibrium spin current in the 
normal region can be interpreted as follows. Electrons of 
the normal 2DEG, while they are scattered back by the 
SOI region, suffer the action of an effective magnetic field 
which forces their spin to precess. This field is directed 
normally to the wave vector of the backscattered wave. 
It is easily seen that due to the precession the spin vec- 
tors of both polarizations of incident electrons get equal 
augmentations of the x projection while the variations 
of other projections are of opposed signs. Electrons in- 
cident with opposed ky produce the spin projections on 
the X axis of opposed sign, so that the total spin density 
is zero. However, the spin currents produced by electrons 
with opposed ky are added. 

In conclusion, the edge spin currents exist in 2D sys- 
tems with spatially nonuniform SOI at thermodynamic 
equilibrium. These currents exist even in the normal 
2DEG if it borders with a SOI region. The spin den- 
sity flows tangentially to the border within a strip of the 
width, which is determined by the characteristic wave 
vector of the SOI, the spins being polarized in the plane 
perpendicular to the border line. These spin currents are 
proportional to the SOI constant a, and are much higher 
than the equilibrium spin current in an infinite 2D SOI 
system, since the latter is proportional to a"^. The edge 
spin currents originate from the spin-dependent scatter- 
ing of electrons from the nonuniform SOI region. 
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